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Abstract 

The classical Bernoulli and baker maps are two simple models of deterministic chaos. On 
the level of ensembles, it has been shown that the time evolution operator for these maps 
admits generalized spectral representations in terms of decaying eigenfunctions. We intro- 
duce the quantum version of the Bernoulli map. We define it as a projection of the quantum 
baker map. We construct a quantum analogue of the generalized spectral representation, 
yielding quantum decaying states represented by density matrices. The quantum decaying 
states develop a quasi-fractal shape limited by the quantum uncertainty. 

1 Introduction 

It is our privilege to contribute a paper to the memory of Prof. Shuichi Tasaki. On a few 
occasions Prof. Tasaki encouraged us to publish the present work. We devote this work to him. 

The present work concerns the subject of quantum chaos. More precisely, the study of 
quantum systems whose corresponding classical systems exhibit chaotic behavior. Prof. Tasaki 
was one of the first authors to link irreversibility and chaotic dynamics using "generalized spectral 
representations" of time evolution operator (also called Frobenius- Perron or FP operator). One 
interesting outcome of the work of Prof. Tasaki and others (Refs. [I]- [5], as well as [6] and 
references therein) was the demonstration that eigenfunctions of the FP operator may have a 
fractal nature. This was shown using classical chaotic maps, such as the multi-Bernoulli map or 
the multi-baker map. The present paper is motivated by this work of Prof. Tasaki and others. 

In classical systems, chaos may appear in isolated systems with few degrees of freedom as 
a consequence of stretching and folding dynamics. A simple model of this type of chaos is the 
baker map. The baker map acts on a unit square with coordinates {q,p) representing the phase 
space. The square is squeezed down in p direction; it is stretched in q direction by a factor of 2 
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Figure 1: The classical baker map. 



and then right half is put on top (see Fig. [T]) . These time-evolution rules are a simple example 
of stretching and folding dynamics. Their consequence is a chaotic time evolution where any 
uncertainty in the initial condition grows exponentially with time until the uncertainty is spread 
over the whole phase space. 

Despite its chaotic evolution, the baker map is invertible and unitary. After applying the 
map any number of times, the resulting final phase-space distribution can be reverted to the 
initial distribution by application of the inverse map. 

An even simpler map is obtained by projecting the baker map onto the horizontal (q) axis. 
This simpler map is called the Bernoulli map; in contrast to the baker map it is not invertible. 
The Bernoulli map maps a number Qt G (0, 1] as 



where r is a discretized time variable, taking only integer values. 

One can understand the chaotic nature of the Bernoulli map by expressing the variable q 
in binary notation. Written this way, the effect of the Bernoulli map is simply a digit shift. If 
Qr = 0.did2d3-.- then qr+i = 0.^2^3(^4. where = or 1. An initial q given by an irrational 
number with random digits will be mapped to a random new point; the trajectory followed by 
q over time will be a random trajectory. 

Corresponding to this behavior of a single trajectory, an ensemble of trajectories initially 
close together will end up uniformly spread out over the whole phase space (i.e., over the 
whole unit interval). An ensemble uniformly distributed over the whole phase space is an 
equilibrium distribution, since further applications of the map will not change it. Equilibrium 



Qr+i = rnod 1. 
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is an attractor for ensembles; any initially smooth distribution of trajectories will approach the 
equilibrium distribution. This situation is reminiscent of the irreversible approach to equilibrium 
of ensembles in macroscopic systems, exhibited in processes such as diffusion. 

In fact, as discussed in Refs. [H El [6] irreversibility (or time-symmetry breaking) can be 
formulated explicitly in simple chaotic models, such as the Bernoulli and baker maps. In precise 
terms, irreversibility means that the FP operator for ensembles admits eigenvalues with mag- 
nitude less than 1, which implies that the corresponding eigenfunctions decay in time. Such 
eigenvalues appear within "generalized spectral representations" of the FP operator. A remark- 
able feature of these representations is that they may include decaying eigenfunctions even in 
cases where the FP operator is unitary, as is the case in the classical baker map. These repre- 
sentations involve both regular functions and generalized functions (or functionals) playing the 
role of right and left eigenstates of the FP operator, respectively. 

In quantum mechanics, chaos is a more involved concept. For small, isolated quantum 
systems the Hamiltonian has a discrete real spectrum. The eigenvalues of the FP operator are 
complex numbers with magnitude equal to 1; no approach to equilibrium is apparent because 
there are no decaying eigenfunctions. 

To describe quantum chaos one can introduce certain criteria, such as integrability, or the 
spacing of energy levels [7] (a thorough discussion on integrability of quantum systems is given 
in Ref. [8]). In this paper, we will follow a different approach. We will consider quantum maps, 
where the time evolution is applied in steps (see Refs. [11^ \T2\ [T3] for a discussion on quantum 
stochastic maps). We will study density matrices, corresponding to the classical statistical 
ensembles. And we will assume the quantum maps are coupled to their environment, which 
causes instantaneous loss of coherence at each time step [U [10] . 

Specifically, we will study the quantum baker map coupled to its environment. The quantum 
baker map is the quantum analogue of the classical baker map. It has a physical realization as 
a set of quantum bits that can be shifted by an external action [T3], one step at a time. We will 
assume that the environment causes the quantum bits to experience instant decoherence after 
each application of the quantum map. The result of the combination of the quantum baker map 
with decoherence will be a new map, which will turn out to be the quantum analogue of the 
classical Bernoulli map. We will call this new map the quantum Bernoulli map. 

We will show that in the quantum Bernoulli map density matrices approach equilibrium, 
similar to their classical counterparts. The approach to equilibrium will be described by in- 
troducing quasi-eigenstates of the quantum FP operator, which decay in time. These states 
will appear in an expansion of the quantum FP operator analogous to the generalized spectral 
representation of the classical FP operator. 

The present work is of interest because it exemplifies the quantum correspondence between 
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classical chaotic maps and their quantum analogues using an approach that as far as we know 
has not been taken before, namely through the generalized spectral representation of the FP 
operator. We will show that the quasi-eigenstates of the quantum FP operator approach their 
classical counterparts in the classical limit. Interestingly, when approaching this limit, the 
eigenfunctions take a self-similar, quasi-fractal shape up to the point allowed by Heisenberg's 
uncertainty principle. (Fractals and self-similarity in quantum systems have been studied in Refs. 
|15j-|20|. although in different contexts. A connection between the decaying eigenstates of the 
classical FP operator and spectral properties of chaotic quantum systems has been investigated 
inRef. [21]). 

In the following section we will summarize the ensemble description of the classical Bernoulli 
map, including the formulation of decaying eigenstates of the FP operator. Then we will describe 
the quantum Bernoulli map and we will formulate the quantum analogue of the generalized 
spectral representation of the FP operator. We will finish by discussing the development of the 
quasi-fractal shape in the quasi-eigenstates and giving final remarks. 

2 The classical Bernoulli map 

In this section we will review the classical Bernoulli map, following Ref. [6]. We will focus on 
an ensemble description of the maps, where the central role is played by the FP operator and 
its spectral (eigenvalue) representation. As we will see the right eigenstates of the FP operator 
are smooth functions that decay in time, while the left-eigenstates are functionals. The set of 
right and left eigenstates together with the corresponding eigenvalues will form the generalized 
spectral representation of the FP operator. 

To introduce the Bernoulli map we start with the baker map illustrated in Fig. [TJ Its FP 
operator is defined as 

P('/A2p), 0<p<l/2 

^(^/2 + l/2,2p-l), l/2<p<l. ^2) 

for a probability density p{q,p) defined in the unit square. The Bernoulli map is a projection of 
the baker map, obtained by integrating over a re-scaled p, 

+ 1 



This gives the FP operator for the Bernoulli map as 

l + q 



UBp{q) = \ 



+ P 



P] (3) 



(4) 



.2/ V 2 

where p{q) is a probability density defined in the unit interval. Successive applications of Ub 
on p{q) will make it evolve towards the uniform equilibrium distribution p°'^{q) = 1. This is 



4 



true provided p{q) is a normalizable function, which excludes delta-functions corresponding to 
trajectories. To define normalizable functions we will use the inner product 

{X\Y)= [\qX*{q)Y{q) (5) 





where Y{q) = {q\Y) and X*{q) = {X\q). The norm of p is then {p\p). 

We are discussing the approach to equilibrium of probability densities. But there is a property 
of the Bernoulli map that seems to indicate that there can be no approach to equilibrium at 
all! Solving this apparent contradiction will lead us to the main result of this section. We first 
introduce the Hermitian conjugate operator J/jj in the usual way 

X\uI\y) = {Y\Ub\XY (6) 



u 



This equation shows that the right eigenstates of f/jj must be left eigenstates of Ub and vice- 
versa. Using the explicit form of Ub the explicit form of can be shown to be 

^B^W-<^ Y{2q-1), l/2<q<l. 

The property of the Bernoulli map that concerns us here is that while C/jj is not unitary because 
U^B^B / 1) it is isometric: 

UbUI = 1 (8) 

This is an important property because it means that any right eigenfunctions of that 
are normalizable must have eigenvalues with magnitude equal to one. To see this, note that 

p\UbU^b\p) ~ ip\p) ~ ^ P with unit norm, including any normalizable eigenstates of 

t 

B- 

The isometry of f/jj implies that it has no normalizable, decaying right-eigenfunctions, be- 
cause normalizable eigenfunctions must have eigenvalues with unit magnitude. This, in turn, 
implies that there are no normalizable, decaying, left eigenfunctions of Ub- All this seems to 
negate the approach to equilibrium characterized by decaying eigenfunctions. However, there is 
a way out. As shown in Refs. [HEllG], Ub can have decaying left eigen- junctionals. In fact, the 
functionals (with a a non-negative number) defined by 

are left eigenstates of Ub with eigenvalues 2~", which means that 

{bJ^Ub = ^(Bo\. (10) 
Taking the Hermitian conjugate we obtain 

^B \Bo) = ^\B^ , (11) 



which means that \Baj are right eigenstates of C/jj with eigenvalue 2 Thus ahowing func- 
tionals to play the role of eigenfunctions allows for the existence of eigenvalues smaller than 
1. 

The operator Ub also has eigenvalues with magnitude less than 1. However, the correspond- 
ing right eigenfunctions are normalizable. Indeed, since U^^Ub 1, the argument below Eq. 
([8]) does not apply to Ub- As a result, Ub can have normalizable right eigenfunctions with 
eigenvalues with magnitude under 1. We will discuss these eigenfunctions now. 

The right eigenfunctions of Ub can be constructed by noting that U b acts as a shift operator 



and 



Uie.M = ] ; (13) 



on the Fourier-basis functions 

ej-/(g) = exp[27ri2J(2Z + l)q] (14) 

with j > and — oo < / < oo integers. 

The existence of shift states allows the construction of a type of coherent eigenstates of Ub, 
given by 

^"(^)"-|J[2vrz2.(2! + l)]"^^''(^) 
It turns out these are the Bernoulli polynomials (q) of degree a [H [U [6] . We have 

Ub\B^) = ^\B^) (16) 

Summarizing, while the right eigenstates of Ub are smooth, normalizable functions (polyno- 
mials), the right eigenstates of U^ (and thus the left eigenstates of Ub) are non-normalizable 
functionals. In order to obtain eigenvalues of C/jj with magnitude less than 1 we had to leave the 
domain of regular functions. This is an interesting point. It ties irreversibility to a mathematical 
formulation involving extended function spaces. 

The Bernoulli polynomials and the functionals ( Ba are orthogonal, 



B^\Ba)=6p,a. (17) 
As a result, in the domain where J2a \^a) ~ operator can be represented as 

oo 

W9) = E2"°"^a(g)(S,|p) (18) 

where BQ{q) = 1. This is the generalized spectral representation of the FP operator. The 
representation (jlSp displays the decay rates 2~", which are powers of the Lyapounov exponent 
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exp(ln2) of the Bernoulli map. This representation clearly shows that when r — >• oo probability 
densities approach the uniform, equilibrium density p*^^ = Bq^q) = 1. As discussed in [H [5l 
[B] the representation (jlSp is valid for smooth, differentiable probability densities (which can 
be expressed as superposition of Bernoulli polynomials). As already mentioned, this excludes 
trajectories, which are represented by Dirac delta functions. 

Setting r = in Eq. (fT8|) and evaluating the inner product (^Ba\p^ explicitly we obtain, 
with = d'^p/dq°', and po = (BqIp) 



oo -I 

Pil) = P0 + Y.-Mq)\p^''-'\i)-P^''-'\0)] (19) 

This is another form of the generalized spectral representation (llSp for r = 0. It is also known as 
the Euler-Maclaurin expansion. It is the generalized spectral representation of the unit operator 
U]j = 1 for T = 0. For later comparison with the quantum Bernoulli map we write this expansion 
as 

piq) = PO (20) 

+ E ^ [i?«(9)p("-i)(i) - (-i)"i?,(i -(/)p("-i)(o) 

a=l 

where we used the symmetry property of the Bernoulli polynomials, 

B^{q) = i-irBM-q)- (21) 
For T > the generalized spectral representation (llSp takes the form 



Uhp{q) = PO (22) 
1 1 
2^ a! 



oo 



+ E \B^{q)p^-'Hl) - {-ITBM - '7)p("-^nO) 



a=l 

3 Quantum baker map 

In the following we will focus our attention on the quantum version of the Bernoulli map and 
the representation corresponding to Eqs. (jlSp or (j22p . Before going to the quantum Bernoulli 
map, in this section we will give a brief introduction to the quantum baker map, upon which 
the quantum Bernoulli map will be built. 

We will follow the definition of the quantum baker map given in Ref. [22j . As in the 
classical case, we have a unit square, with horizontal coordinate g„ ("position") and vertical 
coordinate pm ("momentum"). Since this is a closed, finite quantum system, these coordinates 
are quantized. We divide the unit square into an N x N grid. We assume that N = 2^ where 
D is an integer. The quantized position and momentum are given by 

Tl 

qn = j^, n = 0,l,2...iV-l (23) 
7 



771 

Pm = ^, m = 0,l,2...iV-l (24) 

The position states or the momentum states \pm) are a basis for all possible states of the 
system. They satisfy the orthogonality relations 

{in I Qn') = 5n,n', {Pm \ Pm') = Sm,m' (25) 

and completeness relations 

N-l 



El'?")(9n| = l (26) 
n=0 

E \Pm) (Pml = 1 (27) 
m=0 

The transformation from the position to momentum representation is given by 

1 

{qn I Pm) = —7= eyip[{ipmqn)/K\ 

V N 



^ exp[{27rimn)/N] (28) 



/N 
where 

n = ^— (29) 
27rN ^ ' 

plays the role of Planck's constant in this model. The classical limit is given by — )• oo, 

corresponding to an infinitely fine grid. 

The quantum baker map is represented by a unitary operator T. In position representation 

its matrix elements are 



V2 



Em=o ^ exp [2TTim{n' -2n)/N], < n < N/2 



{qn' \T\qn) = ^x{ (30) 



N 



States IV't) are transformed as 



exp [2Trim{n' - 2n)/N] , N/2<n<N. 



\i;r+i) = T|Vr) (31) 
Hence density operators pr are transformed as 

Pr+l = TprT^ (32) 

This defines the FP operator for the quantum baker map. In the classical limit the quantum 
baker map goes to the classical baker map |22j . 
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4 Quantum Bernoulli map 



We will assume that the quantum baker map experiences instantaneous decoherence in 
representation. As a result, the off-diagonal terms of the density matrix will disappear instanta- 
neously after each time step. Therefore, we will only keep the diagonal terms at each step of time 
evolution. The two-dimensional quantum baker map will be projected into a one-dimensional 
map, which we will call quantum Bernoulli map. 

The FP operator for the quantum Bernoulli map is thus defined as 



Pr+l = U^Pt = [TPtT^^ 



diagonal 

= E E Wn') {Qn'\T\qn) {qn\ Pt \ q-n ) {qn\T^ \qn') {qn'\ (33) 
n=0 n'=0 

The quantum Bernoulli map acts on the diagonal component of density matrices, and transforms 
them into new diagonal density matrices. 

To see the quantum-classical correspondence of the Bernoulli map we start by defining shift 
states as 

N-l 

^hl = E \(ln)(ij,Mn){qn\ (34) 
n=0 

where ej^i{x) is defined in Eq. Due to the discrete nature of the coordinate g„, the functions 

efc(<?n) = ej,i (<?«), k = 2^{2l + l) (35) 

satisfy the relation 

ekiqn) = ek±N{qn)- (36) 
Hence from now on we restrict the domain of k (and thus of j and /) in such a way that 

- N/2 <k< N/2. (37) 

This domain goes to the classical domain in the limit N ^ oo. 

Note that because the number N of modes k is finite, the number of shift states is finite. 
As a result, the quantum Bernoulli map is not isometric. The limited number of shift states 
prevents the construction of exact eigenstates of the FP operator as was done for the classical 
Bernoulli map. However, we will find quasi-eigenstates of the quantum Bernoulli map, which 
will approach exact eigenstates in the classical limit. 

For k = 0, the shift state eo(gn.) = 1 is just a constant and we have 

(U^yeoiqn) = eoiqn). (38) 
For A; / we use Eqs. ([30]) and ([Ml) to obtain [23] 

(U^YeM = I '^-^''^'"^ n^'=° - ^ (39) 

I J 



where 

, I 1, for n even 

= I 1 - 2/+i|2Z + l|/iV, fornodd. ^^^^ 

The result ([39]) says that FP operator acts as a weighted shift on ej^i{qn)- To obtain a non- 
weighted shift we define the new shift states 

ej,i{qn) = fij,iin,N)ej^i{qn), (41) 

where 

Hj^i{n,N) 
The new shift states satisfy 

(uSVhM = { ^f-''""'' ] ; (43) 

Up to this point we have a complete quantum-classical correspondence (compare Eqs. (jl3p and 
()43p ). The weight factor in Eq. (I40p has the form 



1 

j'=i 



(42) 



Sjy(n,iV) = l + 0(;i) (44) 
so it goes to 1 in the classical limit (i.e., N ^ oo with l,j finite) and we recover Eq. (I13p . 



5 The quantum analogue of the generaUzed spectral represen- 
tation for r = 

We will study now the quantum-classical correspondence of the Bernoulli map based on the 
generalized spectral representation discussed in Section [2l We will first find an expansion of the 
unit operator {U^Y with r = 0, analogous to the Euler-Maclaurin expansion in Eq. ()20p . In 
the next section we will consider the case r > 0. 

Consider an arbitrary density p{qn)- The uniform part is given by 

^ N-l 
n'=0 

The non-uniform part, 

h{<ln) = P{qn) - PO (46) 



is expanded in terms of the shift states ej^i{qn') as 

^ iV-l [N] 

Spiln) = J^Y^ Y.^jAQn)eliiqn')p{qn') (47) 



n'=0 j,l 
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where the superscript [N] indicates the restriction in Eq. (j37p . In the classical case the gener- 
alized spectral representation involves derivatives. Seeking a correspondence, in the quantum 
case we introduce differences. We use the notations 



a > 



(48) 



and 



pW(g„,,iV) = p(g„0 (49) 

for the differences. In the classical limit, with the condition a <ti N, the differences go to 
derivatives, 

hm p(-)(q,N) = f^p{q) (50) 

Note that the difference at a point gives the function at that point minus the function at 
the point q-n'-i on the left. We could as well have taken the difference between the points qn'+i 
and q^i. As we will see, our definition of difference breaks the symmetry of the eigenfunctions 
around 1/2 (the mid point between the end points and 1) that exists in the classical model. 
The symmetry is restored in the classical limit. 

The summation over n' in Eq. (j47p may be written as 



N-l _ 
R0=Y1 (^*jMn')p{Qn') = {P^^N-l) - p(i?o)} + ^1 

n'=0 ^ e„.,(c/i) 



1 



where 



N 



-1 N~l 



E ^U<ln')P^^\qn',N) 



The summation in Ri is similar to the left hand side of Eq. ()5ip . Thus we have 



[p^^\qN^l,N) - p^'\q,, N)eli{qi)} + R2 



where 



i?2 



N-1 



l-e*,(gi) n'=2 



E4/(9n')P^'^(9n',iV) 



(51) 

(52) 
(53) 
(54) 



By recursion we get 



N-1 



N 



E eli{qn')p{qn') = E 



_jy-{a-l) 



n'=0 



1 i^-^lMi) 



X [p("-^)(g^_i,iV) - p(°-i)(g„_i,iV)e*,(g«_i) 



(55) 



Hence, the density can be written as (see Eqs. ([36|) and (jI7|) ) 

N [N] 



p{qn) 



1 



p^^-^\qN-uN) - p^^-^\qo,.i,N)elMo.-i) 



(56) 
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Let us now define the quantum Bernoulli polynomials 



-ej,i{qn) 



(57) 



(l-e*,(l/iV)) 



a 



(in Appendix Rl we show that these are polynomials of degree a). With this definition, and 
re- writing the last term of Eq. ([55]) in terms of Sq(1 — Qn+i, ^) we arrive to 



This expression is the quantum version of the classical spectral representation written in the 
form of Eq. (j20p . It is a quantum version of the Euler-Maclaurin expansion. Note that, in 
contrast to the classical expansion, the quantum expansion is not symmetric with respect to the 
exchange of the points x = and x = 1 (corresponding to n = and n = N — 1). This is due 
to the difference operation we have used. It might be possible to use a symmetric difference 
operation, but this will not be investigated here. 

6 The quantum analogue of the generahzed spectral represen- 
tation for r > 

In this section we will obtain an expression describing the time evolution of density matrices 
produced by successive applications of the quantum Bernoulli map. We will obtain an expansion 
analogous to Eq. (I22p . which will involve the time-evolved quantum Bernoulli polynomials. 
Strictly speaking our expansion will not be a spectral representation of the quantum FP operator, 
because the states we will construct are not eigenstates of the FP operator. However, these states 
will approach the classical eigenstates in the classical limit. 

The time-evolution of the quantum Bernoulli polynomials is obtained from Eq. (I43p as 




(58) 



1 r 

+ Sa((/n,iV)p(""')(Q7V-l,iV)-(-l)"i?a(l-9n+l,iV)p("-'^(ga-l,iV) 

, a! L 



Q=l 




[TV] 



(59) 



Changing j — t — t- j and using Eq. (I4ip this becomes 




[TV/2-] 



(60) 



and 



{U^YB^{qn,N)=0 for2^>iV 



(61) 
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Hereafter we consider the case 2"^ < N. Using fij^T-^i{n, N) = fij^i{n, N/2'^) and e*_^_^i{l/N) = 
e* ;(2^/iV) we obtain 

{U^rB^{qn,N) = ^B^{qn,N,T) (62) 

where 

' ^" W2-)" ir /^.v(-,iV/2-) (l-e*,(2Viv)" ^ 



are the time-evolved quantum Bernoulh polynomials. In this expression we have (see Eq. (jlO 

^j^r.l{n,N /2'^) _ \ 1, forneven 



fXj,i{n,N/2^) 1 l + 0{2^k/N), fornodd. 



(64) 



with k = 2^ [21 + 1). Thus for n even the time-evolved Bernoulli polynomials in Eq. (j63p have 
the same form as the initial polynomials in Eq. (|57p . except for the change N — )• N 127 . For n 
odd there is an explicit correction coming from the ratio of weight factors in Eq. (|64p . 
Let us introduce the re-scaled variables 

n' = n/2^, N' = N/2'' (65) 

then with q'^, = n' /N' = qn = n/N we have 

Ba{qn,N,T) = Ba{q'ni,N'), forn cvcn (66) 
= -60,(5^/, A^') + "quantum" correction, fornodd 

where the "quantum" correction comes from the second line of Eq. (j64p . This correction will 
play an important role in the next section when we will discuss the quasi-fractal shape developed 
by the time-evolved quantum Bernoulli polynomials. 

Inserting Eq. (j62p into Eq. (158p we get the following expression for the time evolution of the 
density matrix, expressed in terms of the time-evolved quantum Bernoulli polynomials: 

N 

, 2™ a\ 

a=l 

- (-l)"B„(l-g„+i,7V,r)p("-i)((7,„i,iV)] (67) 

This equation is the quantum analogue of Eq. (I22p . To see how the classical limit is reached, 
we write an alternative form of the quantum Bernoulli polynomials 

, , exp(7rinO 
Ba{qn',N) = (68) 

Ajf /O 

^ ^ exp [2TTik{n' + a/2)/N'] + (-1)° exp [-2TTik{n' + a/2)/N'] 



{U^Ypiqn) = po + J2^^\Bo.{qn,N,T)p("-'\qN-i,N) 



f,^^ [2iN'sm{7:k/N')f 
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We may approximate 

2iiV'sm(7rA;/A^') ^ 2'Kik (69) 

provided that 

N' > 1. (70) 

In this case the Fourier components of the Bernoulh polynomials are dominated by small k 
components with k <^ N' . Moreover let us assume that 

a < iV' (71) 

Under these conditions the Fourier components for small k are independent of N' because 
N' appears only in the expression n'/N', which is equal to n/N. This means that we can 
approximate 

Bo,iq'^,,N')^ Ba{qn,N) (72) 
and therefore, neglecting the correction in Eq. (j66p . we get 

B^{qn,N,T)^Baiqn,N) (73) 

Inserting this into Eq. (j62p we get 

{U^yB^{q^,N) « ±-B^iqn,N) (74) 

which corresponds to Eq. (jl6p . So the quantum Bernoulli polynomials behave as the classical 
ones, i.e., as decaying eigenstates of the FP operator. Furthermore have (with a finite) 

IV' /9 

V o// Ar/N r exp[27rifcn7Af'] + (-l)"exp[-27rffcn7Af'] 

hm B^{q'^,,N') = Inn ^ ' ' ^ 75 

The right-hand side is an expression of the classical Bernoulli polynomial Ba{qn)- 

At r = 0, provided only small values of a with a <ti N contribute to the summation in Eq. 

(jSSp . the discrete difference goes to the continuous derivative in the limit iV — t- oo (see Eq. ()50p ). 

We recover the classical Euler-Maclaurin expansion (j20p . For r > and — )■ oo in Eq. (j67p we 

recover the generalized spectral representation in Eq. ()22p . 

We remark that the existence of quantum decaying states may be interpreted as a signature 

of quantum chaos. The decomposition (|67p shows that any density will approach equilibrium 

with the decay rates 1/2". For 2'^ = N the quantum Bernoulli polynomials vanish identically 

and equilibrium is reached. 
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7 Quasi-fractals 



We consider now the quantum corrections in the Bernouhi map and the development of a quasi- 
fractals shape in the evolving quantum Bernoulli polynomials. 

Let us first now assume that n (in g„) is even such that n' = is an integer. Keeping the 
assumptions ([TDD ^^^'^ GH)) both Ba{qn, N) and Ba{q'^i, N') give a representation of the quantum 
Bernoulli polynomials with different number of grid points, namely N and N' . The point g„ = 
n/N = q'^, = n'/N' belongs to both grids. If both N,N' > 1, then Bo,{qn,N) ?s Baiq'„,,N'). 
After re-scaling by 2*°, the quantum Bernoulli polynomials remain approximately constant at 
the points qn where is integer. Moreover, as discussed in the previous section they are 

approximately equal to their classical counterparts. They behave as decaying eigenstates of the 
FP operator. 

In contrast, if n is odd or more generally, if n' = is not an integer for r > 0, then g^, 

will not belong to the grid with N' points. We expect a deviation from the classical Bernoulli 
polynomial. This deviation is due to the discretization of space, so it will give a quantum 
correction of the order of fi. In addition there will appear the quantum correction in Eq. (j66p . 

Writing n = 2'^(2r + l) with both 7 > and r integers the condition that n' = is integer 
is translated as 7 > t. We denote the set of integers n satisfying this condition as 5,-. We have 

50 = {0,1,2,3,4,5,6,7,8,...} 

51 = {0,2,4,6,8,...} 

52 = {0,4,8,...} 

(76) 

At each time step, the quantum Bernoulli polynomials act as quasi-eigenstates of the FP operator 
only at the points in the sets Sr- At the other points we get deviations. The recursive nature 
of the sets in Eq. (j76p means that these deviations appear in a self-similar fashion. This gives 
the evolving Bernoulli polynomials a quasi-fractal shape. 

In order to visualize this, we have employed a computer program to calculate the exact 
evolution of densities under the quantum Bernoulli map. As an example, at r = we take as 
the initial state the polynomial Ba{qm with a = 3 and = 4096. At each step, we re-scale 
the density by 2'^. If this were the classical Bernoulli map, the graphs would remain unchanged 
for r > 0, because B^{x) is an eigenstate of the FP operator with eigenvalue 1/2'^. But for 
the quantum map we have a different behavior. The result is shown in Figure El We see the 
graphs developing the quasi-fractal shape mentioned above. We also notice that as r increases 
the graphs are shifted. This is due to the a term in the numerator of the right hand side of 
Eq. (|68]) . Provided a <^ N' = N/2'^ the shift is negligible, but as r increases the shift becomes 
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Figure 2: Plot of {2^U^Y B'i{qn, N) for N = 4096 
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Figure 3: Plot of {2^U]^Y B^iqn, N) for N = 4096 and r = 6 
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Figure 4: Zoom-in of Figure [3] 

noticeable. 

Figure [3] shows a zoom-in view of Ba{qn,N) with a = 3, = 4096 and r = 6. Figures 0115] 
show further zoom-in views. 

The self-similarity of the graphs stops when the space resolution is of the order of 1/N ^ % 
(see figure [5]). For this reason, we call this evolving state a quasi- fractal. 

When N — )• oo, the quasi-fractal behaves as a true fractal. At the same time, when — ?• oo, 
the amplitude of the fractal deviations from the classical Bernoulli polynomial gets smaller and 
the graph looks smoother. Eventually, it just looks like the classical Bernoulli polynomial (see 
Figures [6] and C]). 

8 Concluding remarks 

We summarize the main results. We introduced the quantum Bernoulli map by coupling the 
quantum baker map to an external environment that produces instant decoherence after each 
time step. 

We described the quantum Bernoulli map in terms of decaying quasi-eigenstates of the 
Frobenius-Perron operator. These states are analogous to the classical Bernoulli polynomials. 
We found conditions under which these quantum states approach the classical ones. Moreover 
we found a quantum analogue of the generalized spectral representation of the classical Bernoulli 
map. We also found a quantum analogue of the Euler-Maclaurin expansion. In this expansion 
the classical Bernoulli polynomials are replaced by their quantum counterparts; derivatives are 
replaced by differences. 
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Figure 7: Plot of {2^U^Y B-i{qn, N) for N = 4096 and r = 6 (compare with Fig. ED. As N 
increases [fi decreases), the curve looks smoother. At the same time, the self-similar pattern 
appears at smaler length scales. 

We have suggested that a signature of quantum chaos may be the presence of decaying 
eigenstates (or quasi-eigenstates) like the ones we have described in this paper. 

One interesting finding is that even after decoherence, the quantum Bernoulli map shows 
quantum corrections with respect to the classical Bernoulli map. The quantum corrections 
give a quasi-fractal shape to the evolving quantum Bernoulli polynomials. These corrections 
vanish in the classical limit. At the same time, as the classical limit is approached, the quasi- 
fractal approaches a true fractal. In a sense, this fractal is hidden in the classical limit. To our 
knowledge this is a feature chaotic quantum systems that has not received attention before. 

An open question is: do quasi-fractals like the ones we described here appear in other chaotic 
quantum systems? Self-similarity in quantum systems has been discussed in other contexts. For 
example in [TJ [15], self-similarity has been discussed in the context of non-linear resonances. 
Refs. [171 [181 [H] have discussed singular spectra or fractal spectra of quantum systems, that 
appear, for example, with quasi-periodic lattices. Ref. [20] gives an interesting discussion on 
fractals generated by a quantum rotor. The self-similarity discussed in the present work is of 
a different origin. Still, it would be worthwhile to investigate any connections with these other 
works. 

Another question for future work is whether the quantum baker map can be described along 
the lines developed in this paper (i.e. using quantum Bernoulli polynomials and difference 
operators). 

A final point worth noting is that the model we have described may be realized as a set of 
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shifting quantum bits subject to repeated measurements after each shift. The measurements 
cause instantaneous decoherence. If the measurement time interval is greater than the Zeno 
time [21], the system decays exponentially [25]. Equations (j62p and (j74p manifest this prop- 
erty. However if the time interval between measurements approaches zero, we might obtain the 
quantum Zeno effect. We speculate that repeated measurements will prevent the system from 
decaying despite the underlying chaotic dynamics. 
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A The quantum Bernoulli polynomials 

In this appendix we show that the functions 

al -e i{q ) 
Ba,N{qn) = -^Yl T7~V^ C^^) 

are polynomials of degree a in g„. Let us first state two properties of these functions, which are 
easily proved: 

N-l 

E Ba,N{qn) = (78) 



n=0 



and (for n > 0) 



Ba,Niqn) - Ba,N{qn-l) = (a/ N)Ba-l,N (qn) , a>l 

Bi^Niqn) - Bi^N{qn-i) = 1/iV, n>0 (79) 
From the last equation we find that 

Bi,N{qn)=qn + Bi^N{qQ) (80) 
We can find Bi^N^qo) using Eq. (f80]l together with Eq. (|78l) . which gives 

N 

Y,qn + NBi^Niqo) = 0. (81) 

n=l 

Using 

n=l 

we get 

1 1 

Bi,N{qn) = ~ 2 + 2iV ^^^^ 
20 



In a similar way, using 



n=l 



iV3 AT 

~3 2~ "6 



(84) 



we find that 



B2,N{qn) =ql-qn + ^ + -^qn - T7 + 



1 5 



(85) 



We can continue in this way for a = 3, 4, • • • showing that Ba^N{qn) are polynomials of degree 
a. Moreover, recalling that ~ h and consulting a table of classical Bernoulli polynomials 
[6] , we find that 



for a = 1, 2. Similar relations must hold for a > 2 with a <^ N, because of Eq. ()75p . 
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